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Statistical isotropy is often assumed in cosmology and should be tested rigor- 
ously against observational data. We construct simple quadratic estimators to re- 
construct asymmetry in the primordial power spectrum from CMB temperature and 
polarization data and verify their accuracy using simulations with quadrupole power 
asymmetry. We show that the Planck mission, with its millions of signal-dominated 
modes of the temperature anisotropy, should be able to constrain the amplitude of 
any spherical multipole of a scale-invariant quadrupole asymmetry at the 0.01 level 
{2a). Almost independent constraints can be obtained from polarization at the 0.03 
level after four full-sky surveys, providing an important consistency test. If the am- 
plitude of the asymmetry is large enough, constraining its scale-dependence should 
become possible. In scale-free quadrupole models with 1% asymmetry, consistent 
with the current limits from WMAP temperature data (after correction for beam 
asymmetries), Planck should constrain the spectral index q of power-law departures 
from asymmetry to Ag = 0.3. Finally, we show how to constrain models with ax- 
isymmetry in the same framework. For scale-free quadrupole models, Planck should 
constrain the direction of the asymmetry to a la accuracy of about 2 degrees using 
one year of temperature data. 
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I. INTRODUCTION 

In recent years, observations of the cosmic microwave background (CMB) radiation fluc- 
tuations by the Wilkinson Microwave Anisotropy Probe (WMAP) and a large number of 
ground-based and suborbital experiments have led to a precise measurement of the tem- 
perature anisotropy power spectrum up to multipoles of a few thousand ([1-6]) . Apart 
from some claimed "anomalies" (sec below) the observations are consistent with a dark- 
energy-dominated cosmology with statistically isotropic, Gaussian, adiabatic perturbations, 
as expected from simple models of inflation. (We will refer to this as the "concordance" 
ACDM model.) If statistical isotropy apphes, then the harmonic coefficients of the temper- 
ature fleld, 

«L = / dnYC^{n)AT{Q), (1) 

must satisfy 

If the fluctuations arc Gaussian and statistically isotropic, their statistical properties are 
completely described by the power spectrum C^^. 

There have been some hints of "anomahes" in the WMAP data, perhaps suggesting a vi- 
olation of statistical isotropy. These include alignments of low-Z multipoles [7-10], evidence 
for power asymmetry [11, 12] and for a deep cold spot in the southern Galactic hemisphere 
[13, 14]. For the most part, these anomalies have been found by examining the data without 
reference to specific theoretical models. There is, therefore, an a posteriori aspect in comput- 
ing their statistical significance which is difficult to assess [15]. Some authors have, however, 
claimed highly significant discrepancies between the CMB data and the concordance ACDM 
model [16]. 

Interest in the CMB anomahes has motivated theorists to build infiationary models that 
violate rotational invariance, either via the addition of a vector field [17-19] or via an isocur- 
vature perturbation [20, 21]. In addition, a number of phenomenological models that vio- 
late statistical isotropy have been proposed, which can be tested against observations (e.g. 
[17, 22, 23]). 

Reference [17] considers the phenomenological model in which the primordial power spec- 
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trum depends on a preferred direction, 

V{k)^V{k)[l + g{k){k.n)% (3) 

where g{k) is some arbitrary function of wavenumbcr. There has been considerable interest 
in this model recently. The authors of [24] applied Gibbs sampling to the WMAP five-year 
maps to test models with g{k) — g* — constant, finding strong evidence for a preferred 
direction with g^ — 0.15 ± 0.039 using multipoles I < 400. Hanson and Lewis [25] corrected 
some algebraic errors in the analysis of [24] and applied a simpler quadratic estimator to 
the WMAP 5-ycar data. These authors also found evidence for a highly significant (9a) 
departure from statistical isotropy, but with a preferred direction suspiciously close to the 
ecliptic poles, suggestive of some type of systematic effect in the WMAP data. This analysis 
was confirmed by [26] . A subsequent analysis [27] showed that asymmetries in the WMAP 
beams fully account for the observed violation of statistical isotropy. 

Despite this negative conclusion, it is still important to assess the prospects of constrain- 
ing violations of statistical isotropy in the CMB with more precise experiments. In this 
paper, we focus on constraints from the Planck satellite, which was launched successfully in 
May 2009 and has recently completed its third full scan of the sky. The Planck satellite has 
much higher signal-to-noise ratio in both temperature and polarization than WMAP [28]. 
It also has higher angular resolution — 5 arcmin full-width at half-maximum (FWHM) at 
frequencies > 217 GHz — so asymmetries on scales of the beam width should have httle 
effect at multipoles I ^ 1000. It is also expected that the Planck beams will be calibrated to 
high precision from scans of bright planets [29]. Following the successful launch of Planck, 
the European Space Agency (ESA) has approved a mission extension until the on-board 
cryogens are depleted. Planck is therefore expected to produce almost five sky surveys, 
compared to the two sky surveys approved for the nominal mission. This combination of 
high sensitivity, high resolution and extended lifetime allows greater scope for testing for 
systematic effects than was possible with WMAP. For example, it becomes possible to use 
polarization maps independently of temperature maps to test for violations of statistical 
isotropy. 

In this paper, we extend the analysis of [30] to assess how accurately an extended Planck 
mission can be used to test models with an anisotropic primordial power spectrum. This 
paper is organized as follows: In Section II, we summarize some basic properties of the 
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anisotropic model. Section III then applies the quadratic estimator of [25], extended to 
include polarization, to compute forecasts for Planck. Our conclusions are presented in 
Section IV. 

II. THE ANISOTROPIC MODEL 
A. Covariance matrix 

We write the anisotropic primordial power spectrum as 

V{k) = V{k) ( 1 + E 9LM{k)YLM{i^) ) . (4) 
V LM / 

We assume that parity-invariance continues to hold in the mean so that L is restricted to 
even values such that 7'(— k) = Vi^. In this paper, we consider a quadrupole modulation, 
i.e. L — 2, \M\ < 2, with a power-law scale dependence on the wave number QLuik) — 
gLM{ko/ky where the pivot scale ko — 0.002 Mpc~^. Scale-invariant modulation corresponds 
to g = 0. 

The harmonic coefficients of the CMB anisotropy can be expressed as 

j^^Afikmmiik), (5) 

where Af- (/c) are the adiabatic transfer functions, either for temperature {X — T) or E-mode 
polarization (X = E). The primordial curvature perturbation is 7^(k) with statistically- 
homogeneous but anisotropic correlations 

(7^(k)7^*(k')) = (27r)35^(k - k')^P(k), (6) 

with P(k) given by Eq. (4). Thus, the covariance matrix of the harmonic coefficients is 

r^XX' ^ / X X'* 

I\m\,l2m2 \ hmi l2m2 

where 

Qf""' = 47ry" d\nkV{k)Af^{k)Ail'{k) (8) 
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is the usual isotropic power spectrum. The additional term in Eq. (7) due to the power 
asymmetry is 

LM 



x(-l) 



mi 



LM 

(2L+l)(2/i + l) (2/2+1) 
27r 





I L h k) 
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(9) 



where 



Af/f (?) = ^^j d\nkVn{k)^Umf:{k) . (10) 

Since L is even, the anisotropic covariance is nonzero only for li — I2 even as required by 
parity invariance. In the following it will also be convenient to introduce 2x2 matrices 

Ciimi.Zama and Ci^i^ with elements C;^^'^';^^^ and Cfj^''{q) respectively. 



B. Quadratic estimators and the Fisher matrix 

Here we assume that the scale-dependence of the power asymmetry (i.e. q) is known. 
We can then use the quadratic estimator of Ref. [25], extended to polarization, to form 
estimates qlm of the anisotropy parameters. For an isotropic survey and in the limit of 
small primordial anisotropy, these take the form 

9LM — 2 ^LM,L'M' 2^ "Zimi q * "i2m2- y^^) 

L'M' Xiiimi X2hm2 HVM' 

Here, = 'YliX'\i^\°^)^^\'^'^' ^hn ai'^ ^hc temperature and polarization multipoles after 
weighting with the inverse of their isotropic total (signal-plus-noise) covariance matrix. The 
Fisher matrix, evaluated at qlm = 0, is given by 



LM,L'M' 

Z 

hmi l2m2 



(r^tot\-l ^^Iimi,l2m2 ^^tot^-l ^^l2m2,li^i 
OgLM OQL'M' 



(12) 



In the limit of vanishing primordial anisotropy, the inverse of this Fisher matrix equals the 
covariance of the errors on qlm, i-e. -^lm.l/m' = isiLMgl'M')- 

The assumed isotropy of the survey ensures that the Fisher matrix at g^M = is diagonal. 
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Using Eq. (9), the diagonal elements evaluate to 



'LM,LM 



E 

hl2 



Tr(Q,,,(C*f)-^Q,,,(Cr)-') 



A (2/1 + 1) (2/2 + 1) h k L 



Stt 







(13) 



Since qlm is complex for M 7^ 0, we present our simulation results in the next section 

in terms of Glm = V^^gIc/lm) and Gl-m = v^Ini(^LM) for M > 0, and Glo = 9lo- The 
Glm are uncorrelated and have the same variance as the Qlm- 



III. FORECASTS FOR PLANCK 
A. Constraints on the anisotropy amplitude 

In this section, we consider forecasts for the Planck mission. We use the parameters for 
the 143 GHz channel (the most sensitive of the Planck frequency channels) as given in [28]. 
We therefore assume a Gaussian beam with FWHM of 7.1 arcmin and assume uncorrelated 
isotropic noise in the temperature and polarization maps with root-mcan-squarc noise levels 
of cr-r and ap respectively. For the nominal two-sky-survey mission (one year of observation) 
we adopt cry = 12.2 /iK and ap — 23.3 /iK in 3.4 arcmin pixels (Healpix [31] resolution 
-^side — 1024), corresponding to 42 //K- arcmin noise in temperature and 80//K-arcmin in 
Stokes Q and U polarization. We also consider an extended mission of four complete sky 
surveys (two years of observation) with ut and ap reduced by a factor of -\/2- 

To give a feel for the nature of the anisotropy signal, Figs. 1 and 2 show simulated 
sky maps for a noise-free realisation of a scale-invariant {q = 0) quadrupolc-modulation 
model with 5120 = 0.1 {g2M — 0, M ^ 0). These maps have been generated using the 
prescription described in [25], generalized to polarization. This uses an approximate square 
root of the anisotropic covariance matrix, linear in the giM, to simulate maps as the sum of 
a statistically isotropic part and an anisotropic part. The isotropic component of a noise- 
free temperature map is shown on the left in Fig. 1 and the anisotropic component, clearly 
showing a preferred direction along the polar axis, is shown on the right. 

The anisotropic contribution to the Stokes Q and U polarization maps is shown in Fig. 2. 
Here, we have smoothed the polarization maps with a Gaussian of FWHM 3 deg. to en- 
hance the visual impact of the statistical anisotropy in the Q Stokes map. Because of the 
quadrupole asymmetry in the primordial power, modes of the primordial perturbation with 



FIG. 1: Noise- free simulation of a model with scale- invariant quadrupole asymmetry in the pri- 
mordial power with giM = 0.1(5l2'5a/o- The isotropic component of the temperature map is shown 
on the left, and the anisotropic component on the right. The colour scales are in ^K. 
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FIG. 2: Anisotropic components of the Q (left) and U (right) polarization maps in a noise-free 
simulation of the model in Fig. 1. The maps have been smoothed here with a Gaussian beam of 
FWHM 3° to enhance the imprint of the preferred axis in the Q map. The colour scales are in ^K. 

wavevectors along the polar axis tend to have their amplitude enhanced. For such modes, 
the polarization generated by Thomson scattering is pure Q in the polar basis and varies 
in amplitude with the polar angle as sin^ ^ [32]. The dominant effect of the statistical 
anisotropy is therefore observed in the Q Stokes parameter and is concentrated toward the 
equatorial plane. 

To illustrate the machinery summarised in Section JIB, we have generated five simu- 
lations of the scale- invariant quadrupole- modulation model with g2M = 0.16qm and added 
instrumental noise appropriate to one year of observation with the Planck 143 GHz channel. 
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FIG. 3: Estimates of the G20 (left) and (right) anisotropy parameters (shown with points) 

and their (one-sigma) Fisher errors ([Red] sohd hnes) as a function of Zmax from five simulations 
of the model in Fig. 1 for one year of Planck data. The input parameters G20 = 0.1 and G2-2 
are shown with horizontal [Blue] solid lines. From top to bottom we analyse temperature only, 
E-mode polarization only and temperature plus polarization. 

We then estimate g2M via Eq. (11) as a function of /max, the maximum multipole retained 
in the analysis. Since the only nonzero coefficient in these simulations is (720 5 and the survey 
is assumed isotropic, the recovered estimates g2M are statistically equivalent for M 7^ and 
so we show results only for the two (real) components G20 and G2-2 in Fig. 3. We anal- 
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yse temperature alone (top panels), £^-mode polarization alone (middle), and both jointly 
(bottom) . 

With temperature alone, the errors on G2M decrease approximately as l//max over the 
range plotted in Fig. 3 (/ < 1500) reaching 0.005 by /max (in agreement with the minimum- 
variance estimators of [30]). This behaviour follows from simple mode-counting since the 
temperature maps are signal-dominated over this multipole range. ^ However, the polariza- 
tion maps are noise-dominated over much of this multipole range and so the errors approach 
constant values for Zmax ~ 600. Nevertheless, the Planck polarization maps alone can provide 
(almost) independent constraints on an anisotropic modulation to the temperature maps. 
For two sky surveys, the errors on the g2M from polarization arc four times worse than in 
temperature. Consistency between temperature and polarization constraints would provide 
an important test of systematic effects should Planck show any evidence of an anisotropic 
power spectrum. 

B. Constraints on scale-dependence 

In Fig. 4, we compare the Fisher errors on the amplitude of the modulation for a scale- 
invariant model (g = 0) and two models with scale-dependence (g = 1 and q = 2). For larger 
g, the asymmetry in the variance of the Fourier modes is confined to larger scales and so 
relatively more of the constraining power derives from low-Z multipole moments. The low-Z 
modes of polarization are enhanced by scattering at reionization [33] and are expected to 
be signal-dominated in the one-year Planck data. The polarization constraints on the g2M 
therefore become more comparable to those from the temperature as q increases and the 
improvement from observing for longer in polarization lessens. 

If the amplitude of any primordial power asymmetry is high enough, it might be possible 
to constrain the scale-dependence of the asymmetry with Planck. To forecast constraints 
on the spectral index of the power asymmetry, g, we extend the Fisher matrix analysis of 
Sec. II B to include g as a parameter. We must now evaluate the Fisher matrix at nonzero 
Qlm but we assume that the asymmetry is still small enough that we can neglect asymmetry 

^ For modes that are signal-dominated, the scale dependence of the trace term in the Fisher matrix, Eq. (13), 
is weak. Treating the trace as constant gives Fisher information varying as /^g^j which is proportional to 
the number of modes retained in the analysis. 
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FIG. 4: Fisher errors for G20 from temperature, E-mode polarization, and temperature plus po- 
larization in models with power-asymmetry spectral indices q = (left), q = 1 (middle) and q = 2 
(right). For q = and q = 1 we show results for one and two years of observations; for q = 2 we 
show only the one-year errors since they improve very little with further observing time. Note that 
the one- and two-year errors from temperature alone are indistinguishable when q > 0. 



in the (C*°*)-^ terms. The F, 
elements are 



'LM,q 



hmi l2m2 



E 



LM,L'Ar is then unchanged from Eq. (13) but the additional 

1 9^l2m2,hmi 



Ti'(c,,,,(c;f)-'a,Cy,(c;°')-')9LM 



5,9 



hmi l2m2 



f(-^tot\-l dClimi,l2m2 

^^^^^ dq 



E 

hh 



Tr(9,Q,jqf)-^9,Cy,(C,^ 



-itotN-l 



dq 

(2/i + i)(2/2 + i; 

Stt 

1 dCl2m2,hini 

dq 

(2/i + l) (2/2 + 1) 
Stt 




2-1 



E 



X > \9lm\ 

LM 
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(15) 



Note that these elements vanish if qlm = 0. From the scaling of the Fisher matrix elements 
with gLM, we expect the error on q to scale as the inverse of the amphtude of the asymmetry. 

As an example, we consider a fiducial model with scale-invariant quadrupole asymme- 
try with G2M = 0.035a/o- Such a model is compatible with the constraints on quadrupole 
asymmetry from the beam-corrected analysis of WMAP data in Ref. [27], but the ampli- 
tude should be detectable with Planck at the 6a level (fixing g = 0). Forecasts for the 
(marginalised) errors on q for this model are given in Table L With temperature and po- 
larization, Planck should constrain the spectral index to a la accuracy of Aq ~ 0.3. The 
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Aq {la) 


T only 


E only 


T and E 


One year 


0.399 


1.300 


0.322 


Two years 


0.389 


0.878 


0.299 



TABLE I: Fisher errors on the scale-dependence of the power asymmetry assuming a scale-invariant 
quadrupole asymmetry with g2M = 0.03(5mo- 

marginalised errors on the qlm are similar to the case in which q is fixed. 

C. Axisymmetric models 

The preceding analysis makes no assumptions about axisymmetry of the primordial power 
asymmetry. However, if we have good reason to expect axisymmetry, so the model is de- 
scribed by a preferred axis ih and the M — multipoles {q^l} in a frame with the polar 
axis along ih, we can constrain these parameters by post-processing our estimates Qlm- We 
illustrate how this works, assuming a fixed scale dependence for the primordial asymmetry. 

For models with a nearly scale-invariant anisotropy spectrum, many small-scale modes 
contribute to the giM, so we might expect the statistics of the qlai to be approximately Gaus- 
sian. Expressing m in terms of its azimuthal angle a and polar angle ^, m = D{a, /3, 0)z (i.e. 
a rotation of the z direction through Euler angles a and (3), in the Gaussian approximation 
we can write Pr({^LM}|iii, {g*L}) oc exp(— where 

= ^ ^ idLM ~ 9LM)FLM,L'M'{gL'M' " QUM')- (16) 
LM L'M' 

Here, 

qlm = -Dmo(«, 0)9*L (17) 

are the multipoles of the primordial asymmetry rotated from their preferred frame. 
i^MM'i^-i P-il) the Wigner rotation matrices.) If we now assign a uniform prior on 
the direction m [so that ^^{a, ^)dad^ — {A7r)~^ dad cos /3 — (47r)~^dm] and a fiat prior on 
the g^L, Bayes' theorem gives for the posterior 

Pr(m,{g^L}\{gLM})dni = Pr(a, /3, {g^L}\{gLM})dad/3 

oc e-^'/^dui. (18) 
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For an isotropic survey (and assumed weak anisotropy), the Fisher matrix is isotropic 
and we can write Flm,l'm' — Sll'Smm'/ctI- Substituting Eq. (17) into Eq. (16) and using 
the addition theorem for the rotation matrices, we have 



i9lL-29*LY,^[gLMD^Moi(^,/3,0)]+Y,\9LMn • (19) 

L ^ \ M M ) 



Noting that £'mo("> 0) = ^^^^/(SL + l)yLM(m), if we define 

gL (m) = giMYLM (m) , (20) 

M 

(i.e. a map of the estimated qlm at multipole L), we can write the posterior as 

2n 



Pr(m, {g*L}\{gLM})dm oc TTexp - ^ [9*lYlo{z) - hi'^)] 



>^ e^P ( TTTT—TV^ali'i^) ] dm. (21) 



(2L + l)a 



The marginal distribution for the direction of the axis is given by integrating over q^l: 

Pr{rn\{gLM})dm. oc JJexp f (2L+\)a^ ^L(m)^ dm, (22) 

so that contours of constant density are given by the contours of J2l qK'"^) /[(2-^+1)cl]- For 
the simple case of a primordial power asymmetry at a single multipole L, e.g. a quadrupole 
asymmetry, lnPr(m|{^Ljvf}) is proportional to the square of the map of the reconstructed 
multipoles. 

To illustrate these ideas, in Fig. 5 we plot the marginal distributions for the direction and 
amplitude of the model analysed in Sec. Ill A using one of the simulations shown in Fig. 3. 
We parameterise the direction by the Cartesian components of the equal-area projection x — 
2 sin(^/2) cos(q;) and y — 2 sin(/3/2) sin(Q;). With one year of simulated Planck temperature 
data, the constraints on the direction of the axis are Ax ^ Ay = 0.032 (i.e. 1.8°) at 
68% confidence and the x and y components are almost uncorrclatcd as expected from 
isotropy. The amplitude gj,2 = 0.1 ± 0.005 (68% confidence) with the error being very 
nearly (72. With polarization alone, these constraints weaken to Ax w Ay — 0.10 (i.e. 5.7°) 
and 51*2 = 0.11 ± 0.02. The maximum of the posterior in this simulation is at x = 0.042, 
y = -0.035 and g^2 = 0.099 (T only) and x = -0.033, y = 0.062 and g^2 = 0.11 {E only) 
and the at these values are 0.223 and 0.236 respectively, with 5 — 3 = 2 degrees of freedom. 
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FIG. 5: Marginal distributions for the direction (left) and amplitude (5*2, right) from a simulation 
of the nominal (one-year) Planck survey for a model with an axisymmetric quadrupole asymmetry 
aligned with the polar axis {gLM = 0.16l2Smo)- We parameterise the direction with the equal- 
area projection x = 2 sin(/3/2) cos(q!) and y = 2 sin(/3/2) sin(a) and show in solid lines the 68% 
[Red], 95% middle [Blue] and 99% outer [Green] contours from the temperature alone; dashed-line 
contours are from the E-mode polarization alone. For the amplitude (right), we plot the marginal 
distributions from T alone (solid line) and E alone (dashed line). 

IV. CONCLUSIONS 

As summarized in Sec. I, cosmological perturbations are usually assumed to satisfy sta- 
tistical isotropy, as expected in simple models of inflation. Yet there is some tentative ob- 
servational evidence from the CMB suggesting possible departures from statistical isotropy. 
Here we have shown that Planck should be able to set strong constraints on small depar- 
tures from statistical anisotropy and, under favourable circumstance, could set constraints 
on their scale dependence and any preferred direction. 

In this paper we have developed and applied quadratic estimators to test for an asymme- 
try in the primordial power spectrum from temperature and polarization measurements of 
the CMB. Our estimators are optimal in the limit of isotropic primordial power. We have 
tested our methods against simulations that include a quadrupole power asymmetry. 

We have analysed the ability of the Planck mission to constrain models with quadrupole 
power asymmetry using temperature and polarization data. Using temperature data alone, 
Planck should be able to constrain each multipole g2M of a scale-invariant quadrupole 
anisotropy at the 0.01 level {2a), well below the current constraints derived from WMAP 
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(|fl'2M| < 0.07 from Ref. [27]). Using polarization data alone from an extended Planck mis- 
sion (four sky surveys) such an anisotropy can be constrained to an accuracy only about 
three times worse than from the temperature. This offers the possibility of a consistency 
check on the existence on any observed departure from statistical isotropy. 

If the amplitude of a power asymmetry is large enough, it may be possible to constrain its 
scale dependence. We have estimated the Fisher errors when additionally constraining a free 
spectral index describing a power law modulation, QLuik) — gLMiko/ky. For a scale-free 
quadrupole modulation with an amplitude of 1% (i.e. g2o/V^ — 0.01 in an axi-symmetric 
model), we find that an extended Planck mission can constrain the spectral index to a Icr 
accuracy of Aq ~ 0.3. 

Finally, we have considered the constraints on a preferred direction in models with a 
purely axisymmetric modulation of the primordial power spectrum. In a scale-free model 
with a 1% quadrupole modulation, the direction of the preferred axis can be determined 
from Planck data to a precision of about 2° using temperature observations alone and to 
about 6° using polarization data alone. 

The quadratic estimators developed here for isotropic surveys can, in principle, be 
straightforwardly extended to deal with real-world effects such as anisotropic noise and 
Galactic masks. However, the calculation of the in Eq. (11) requires inverse weighting 
the temperature and polarization data with the full covariance matrix for the anisotropic 
survey. This has been done for the WMAP temperature data at its native resolution [25, 34] 
but extending this to polarization and the resolution required for Planck requires further 
work. In practice, fast estimators can still be constructed for anisotropic surveys, with 
moderate loss of performance, by replacing with some heuristically-weighted pseudo 
multipoles (i.e. those computed directly on the masked sky) following techniques used for 
CMB power spectrum estimation (e.g. see [35] for a review). In either case, care must be 
taken to subtract the mean-field response - that obtained on average for no primordial power 
asymmetry - from the quadratic estimator since this is no longer confined to the L — mode 
for an anisotropic survey. The mean-field and the estimator normalisation are then generally 
best determined by Monte-Carlo simulations. 
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